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Abstract 

The basic models of modern cosmology work in FRW spacetime. Hence follows that it is 
important to study the physical and mathematical nature of FRW cosmological models. In 
this work, we consider FRW models with the differential equations of state. Some new classes 
integrable and nonintegrable FRW cosmological models were constructed. It is remarkable 
that all proposed integrable and nonintegrable FRW models admit exact solutions. For some 
of them, such exact solutions are presented. Some artificial two-dimensional FRW models 
were also proposed. Finally in Appendix, we extend the obtained results for g-essence models 
and for its two reductions: k-essence and f-essence. 
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1 Introduction 

FRW models play a central role in modern cosmology. In particular, almost all popular theoretical 
models of dark energy work in FRW spacetime. Again almost all models of dark energy meet 
some difficulties like cosmological constant problems, fine-tuning problems and so on. One of 
consequences of such difficulties of modern cosmology is the necessity more carefully investigate the 
basics of General Relativity (GR), in particular, FRW cosmology. One of the poorly studied sector 
of modern theoretical cosmology is the problem of integrability of cosmological models, especially, 
FRW cosmological models. In our previous papers [1]- [3], we have studied the relationship between 
the basic cosmological equations, namely, the Friedmann equations 

p = -2H-3H 2 , (1.1) 
p = 3H 2 (1.2) 

and some important equations known from the other branches of physics and/or mathematics. 
These last equations include: 
a) The Ramanujan equation 

P = d.3) 
E = (1.4, 

J = ^ (1-5) 



b) The Chazy-III equation 

c) The Lorenz oscillator equation 





v = 2yy - 


3y 2 . 




= o{Y 




Y N 


= X(5- 


-Z)-Y, 


Zn 


= XY - 


- pz. 



(1.6) 

(1.7) 
(1.8) 
(1.9) 
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d) Painleve equations. These six equations, traditionally called Painleve I- VI, are as follows: 



Pj — equation 


y = Qy 2 + t 


Pi I ~ equation 


y = 2y 3 + ty + a 


Pin — equation 


y = j i y 2 -Hy- a y 2 ~f 3 ) + "?y 3 + l 


Piv ~ equation 


y=±y 2 + 1.5y 3 + 4ty 2 + 2{t 2 -a)y+f i 


Py — equation 


y = (^ + ^ 2 -i(y-7y) + ^(^ + f) + ^i) 


Pyi — equation 


y = am + + p + a*)} 



(1.10) 



Here ^(z) = O^ 1 + (y - l)" 1 + (y - z)" 1 ], = [z" 1 + (z - l)" 1 + (y - z)" 1 ], C(*) = 

[7(z-l)(y-l)- 2 + <5z(z-l)(y-z)- 2 ]. 

In this work, we study the FRW cosmological models with the differential equations of state 
(EoS). To construct integrable and nonintegrable reductions of such FRW cosmological models we 
use the implantation method. In particular, to construct integrable FRW models in one dimensions 
we use six Painleve equations (see also pQ- [3])- 

The paper is organized as follows. In section 2, we present the basic equations of FRW cos- 
mology. Section 3 is devoted to study the FRW models with the half EoF and Section 4 to the 
FRW models with the full EoS. In the next section 5, we present the so-called D - models. Some 
artificial two-dimensional models were constructed in the section 6. The last section 7 is devoted 
to the conclusion. 

2 Basic gravitational equations 

We start from the classical GR case. In this GR case, the standard gravitational action has the 
form 

(2.1) 



S = J V^d^xiR + L m - A), 



where R is the scalar curvature and L m is the Lagrangian of the matter. We work with the FRW 
spacetime which has the metric 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



1 - kr 2 



r 2 dfl 2 



(2.2) 



where a(t) is the scale factor, k = —1, 0, 1 represent the three-dimensional space with the negative, 
zero, and positive spatial curvature, respectively. In this case the Ricci scalar reads as 



R = 6 



k 



(2.3) 



where a dot represents differentiation with respect to t. The Friedmann equations read as 



A k 8ttG 



A AttG 



ft- 



3 a 3 



-(/> + 3p), 



(2.4) 



where we note that p < —p/3 implies repulsive gravitation if A = 0. Recaling the Hubble parameter 
H = da -1 these equations can be written as 



rr9 8?rG k A . , 

H 2 = — P-— + -, H — —AitG(p + p) 
3 



k 



3 a^ 



(2.5) 



3 



If the FRW spacetime is filled with a fluid of energy density p and pressure p, then the conservation 
law coud be derived from the Friedmann equations (2.5) as 



p + 3H(p + p) = 0. (2.6) 

In this work, we consider the case: k = A = and set 8ttG = 1. So finally the equations for the 
action (2.1) we can write in the H-form 

p = -2H-3H 2 , (2.7) 

p = 3H 2 , (2.8) 

p = -<3H(p + p), (2.9) 



in the N-form 



p = -27V -37V 2 , (2.10) 
p = 3JV 2 , (2.11) 
p = -3N(p + p) (2.12) 



or in the a-form 



P = (2-13) 
P = 3^J, (2.14) 
p = ~(P+P), (2-15) 

where N — In a. Note that one equation among three is redundant. Here we may choose the first 
two Friedmann equations as two relevant equations. To find H we note that the general solution 
of the equation (2.7) we can write as 



H(a) = y-«~ 3 / Pa 2 da. (2.](>) 

This formula tells us that the EoS has the form 

p = -3a" 3 J pa 2 da (2.17) 

or 

P=-p-3- 1 a Pa (2.18) 

so that its parameter takes the form 

w = -l-3" 1 a(lnp) a . (2.19) 



3 Differential equations of state: 0.5 - models 
3.1 Integrable 0.5 - models 

In this subsection we consider some integrable 0.5 - models that means models with the half EoS 
(half = 0.5). As the integrable cells (which we inplant to the body of the original gravitational 
system of equations), we use Painleve equations. The function of these integrable cells are to 
convert the original gravitational (nonintegrable) system into the integrable system. 
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3.1.1 p - models 

Let's go further or, if exactly, go back. To the Friedmann equations e.g. (2.7)-(2.9) and consider 
the following its extension 

p = -2H-3H 2 , (3.1) 

p = 3iJ 2 , (3.2) 

p = p- 1 p 2 -r 1 (p~ap 2 -P)+ 1 p :i + 5p-\ (3.3) 

p = -3H(p + p), (3.4) 

where a, /?, 7, S are consts. It is the Am a - model (see below). We guess that this system is 
integrable due to integrability of the equation (3.3) which is nothing but the P777 - equation and 
which plays the role of the integrable cell (see e.g. [2]-[l]). To solve it, we start from the equation 
(3.3). In particular, it has the following well-known particular solutions [I]: 

p{t) ee p{t; /1, — pn 2 , A, — Aft 4 ) = ft, (3-5) 
p(t) = p(t; 0, — /1, 0, /ift) = nt, (3-6) 

p{t) = p(t; 2ft + 3, -2ft + 1,1,-1)= t + K (3.7) 

t + k + 1 

p(t) = p(t;p,0,0,-pn 3 ) = K\ft, (3.8) 

p(t) ee p(t; 0, —2ft, 0, 4ft/i - A 2 ) = t[n(lnt) 2 + Xlnt + p], (3.9) 

p(t) ee p(t;-^A,0^ 2 (A 2 -4ft M ),0)= . 2 * , (3.10) 

Kt zv + \t v + p 

p(t) ee p(<;O.5e) = -ei(ln0) f (3.11) 

and so on. Here 

0(t) = f[Ci J„(C) + C 3 Y V (Q], (3.12) 

where = consts, £ = ^/e^t, ^ = 0.5aei and J u (C)-i ^(C) are Bessel functions. We are here 
modest and work just with the simplest solutions. 

i) First let's consider the solution (3.5): 

p(t) = k = const. (3.13) 

From (3.2) we get that 

H(t) =H = 3-°- 5 k - 5 = const, (3.14) 

so that this solution corresponds to the de Sitter case. 

ii) Next, let's consider the other simplest solution namely the solution (3.6): 

p{t) = nt, (3.15) 
where we assume that k > 0. Then the equation (3.2) gives 



i/ = ±^/-i Ub (3.16) 

so that we have 

a = a e ± V¥* 1 ' 5 . (3.17) 
For this particular solution the EoF has the form 

p = -pT3-°- 5 ft a7 V a25 - (3.18) 



The corresponding EoF parameter reads as 

-0.5-0.5.-1.5 



(3.19) 



For this case we have 

n = -\-n.n * I 

3 



a^±a J^t°- 5 e ± ^ t1 - 5 (3.20) 
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and 



3 l V4k j 
Let's calculate the deceleration parameter. We have 



(3.21) 



(3.22) 



Hence we see that q = as to = \f^- ^ means that this solution describes the deceleration 

and acceleration phases of the expansion of the universe and to = \f^i 1S the transition point. 
Similarly, we can explore the other solutions of the Ajjja - model (3.1) - (3.4) as well as the other 
p - models. Now we present the list of p - models. 
1) Ai - models. 



Aja — model 


p = 6p 2 + t 


Aib — model 


Paa =6p 2 +a 


Aic — model 


PNN = 6p 2 + N 


Aid — model 


PHH =6p 2 + H 



2) An - models. 



An a — model 


p = 2p 3 +tp + a 


Ajib — model 


Paa = 2p 3 + ap + a 


Auc — model 


Pnn = 2p 3 + Np + a 


Aiid — model 


Phh = 2p 3 + Hp + a\ 



3) Am - models. 



Am a — model 


p = p- x p 2 - t~ x (p - ap 2 - (3) + jp 3 + dp- 1 


Ahib — model 


p aa = P~ X p\ - a~ 1 {p a - ap 2 - /3)+~fp 3 + SP- 1 


Amc — model 


Pnn = P~ 1 P% - N~ 1 (pN - ap 2 - /3) + jp 3 + Sp- 1 


Ami) — model 


Phh = P~ 1 P 2 h - H~ 1 ipH - ap 2 - /3) + 7P 3 + Sp' 1 



(3.23) 



(3.24) 



(3.25) 



4) Aiv - models. 



Ajva — model 


p = O.Sp- 1 ^ 2 + 1.5p 3 + Atp 2 + 2(t 2 - a)p + Pp- 1 


Aiv b — model 


p aa = 0.5/9- Va + + W + 2(a 2 - a) p + f3 p' 1 


Aivc — model 


Pnn = 0.5p"Vw + 1 V + ANp 2 + 2(N 2 - a)p + pp- 1 


Aivd — model 


Phh = 0.5/9-Vh + 1 V + ^Hp 2 + 2{H 2 - a)p + fip- 1 



(3.26) 



G 



5) Ay - models. 



Ay a — model 


'p = 4>p 2 - t _1 (p - 7p) + t- 2 (p - l) 2 {ap + fip- 1 ) + 5iP 


Ay b — model 


Paa = <t>pl ~ a _1 (Pa - 1 p) + OT 2 {p - l) 2 (ap + fip- 1 ) + 5^ 


Aye — model 


Pnn = <t>p 2 N - N-\p N - jp) + N- 2 (p - l) 2 (ap + pp- 1 ) + Sip 


Ay d — model 


Phh = <t>p 2 H ~ H-\p H - 7 p) + H- 2 (p - l) 2 (ap + pp- 1 ) + 5<P | 



(3.27) 



Here = O.Sp- 1 + (p - V = p{p + \){p - l)" 1 . 

6) Ay i - models. 



Ay i a — model 


p = ^(t)p 2 mp ~ "X^Ha +%+ at)} 


Ay ib — model 


Paa = V {a)p 2 a Z(a)Pa PiP a ^ a % a) + C(«)] 


Ayic — model 


p NN = ^(N) P 2 N £(N)p N - "fcftl-p C(N)} 


Ay id — model 


p HH = V {H)p 2 H aH) PH [a + + C{H)] 



(3.28) 



Here ip[z) = 0.5b" 1 + (p - l)" 1 + (p - z)-\ £(z) = [z' 1 + (z - l)" 1 + (p - z)~% ((z) = 
[ 7 (z-l)(p-l)- 2 + 5z(z-l)(p-z)- 2 ]. 



3.1.2 p - models 

Now we want consider one of important part of integrable 0.5 models namely p - models. Let's 
implant the integrable cell, e.g. the Bins - model (see below), to the body of the gravitional 
equations. The aim of this implantation is the original (in general nonintegrable) system convert 
into the integrable system. If this integrable cell is, e.g. the Bhib - model, as result we get the 
following closed system of the equations 

p = -2H-3H 2 , (3.29) 

p = 3H 2 , (3.30) 

Paa = p^pl-a^ipa-ap 2 - f3)+ W 3 + 5p~\ (3.31) 

p = -3H(p + p), (3.32) 

where a,f3,-y,5 are consts. This system we also call the Bins - model, which is (as we expect 
and believe) integrable. Also we note that the equation (3.31) is the Pm - equation and plays 
the role of the integrable cell. As an integrable, the system (3.29)-(3.32) admits (may be infinite 
number) exact solutions. The construction its exact solutions we start from the equation (3.31). 
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In particular, the equation (3.31) has the following particular solutions (see e.g. [1]): 

p(a) = p(a; fx, — /ik 2 , A, — Ak 4 ) = k, 
p(a) = p(a;0,—n,Q, [in) = na, 

a + k 



(3.33) 
(3.34) 

(3.35) 

(3.36) 
(3.37) 

(3.38) 

(3.39) 

cf>(a) = a v [CiJ v {Q + C 2 Y u (0], (3.40) 

where Ci — consts, £ = y/e^a, v — 0.5aei and J V {£),Y V {£) are Bessel functions. Now we 
consider some simplest solutions. 

i) Let's start from the solution (3.33) 



p{a) 

p{a) 
Pifl) 
p{a) 
p{a) 



= p(a; 2k + 3, -2k +1,1,-1) = 

a + k + 1 

= p(a; fi, 0, 0, — fin 3 ) — Kyfa, 

= p(a; 0, —2k, 0, 4np — A 2 ) = a[«;(ln a) 2 + A In a + p], 

„v-l 

= p(a;-v 2 \,0,v 2 (\ 2 -Akp),0) = 



Ka 2v + \a" + [i ' 



p(a;0.5e) = -ei(ln0) a 



and so on. Here 



p(a) = k = const. 
To find H we use the formula (2.16). As result we get 



H(a) = yZ-a-zfi-iKcP + C) = y/ -Z' 1 k - Ca~ 3 = ^/JF 1 



p a- 



3-!A. 



(3.41) 
(3.42) 



where po — — 3C, A = —k. It is nothing but ACDM model. So this solution corresponds to the 
ACDM cosmology. In this case, the above formulas give 

p = -A, p = p a~ 3 + A (3.43) 

which corresponds the EoF parameter 

w= -l + po^o+Aa 3 )- 1 . (3.44) 

ii) Now let's we consider the next simplest solution namely the solution (3.34): 

p{a) — Ka. (3.45) 

Eq. (2.16) gives 

H{a) = y/-a,- 3 (0.25Ka 4 + C). (3.46) 
The corresponding energy density is given by 

p = -3a- 3 (0.25Ka 4 + C). (3.47) 
For this particular solution the EoF and its parameter take the form 

p = -0.75p - 3CkV 3 (3-48) 

and 

4 

Ka 

w = -ojwT3C' (3 - 49) 

respectively. In the limit a — ¥ oo we get u> — > —4/3. 

Similarly we can write Bj - models induced by the other Pj - equations. Here the list of such 
models. 

1) B/ - models. 



Bia — model 


p = 6p 2 + t 


Bib — model 


Paa =6p 2 +a 


Bjc — model 


Pnn — 6p 2 + N 


Bid — model 


Phh = Gp 2 + H 



(3.50) 
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2) B// - models. 



Biia — model 


p = 2p 3 + tp + a 


Bub — model 


Paa = 2p 3 + ap + a 


Buc — model 


p NN = 2p 3 + Np + a 


Bu d — model 


Phh = %P 3 + Hp + a 



(3.51) 



3) Bui - models. 



Bui a — model 


P = P 1 P 2 — t 1 (p~ ap 2 - f$) + 7P 3 + Sp 1 


Buib — model 


Paa = P _1 Po - a- x (p a - ap 2 - P) + 7P 3 + Sp^ 1 


Bine — model 


Pnn = P~ X P N - N~ 1 (p N - ap 2 - /3) + 7p 3 + Sp^ 1 


Buid — model 


Phh = P^Ph - H~ 1 {p H - ap 2 - ft) + 7p 3 + Sp' 1 



(3.52) 



4) Biv - models. 



Biva — model 


p = 0.5p~ 1 p 2 + l.bp 3 + Atp 2 + 2(t 2 - a)p + ftp' 1 


Bivb — model 


p aa = 0.5j5 _1 Pa + 1 - 5 P 3 + 4a P 2 + 2 ( ft2 ~ a )P + PP~^ 


Bivc — model 


p NN = O.hp^PN + 1.5p 3 + ANp 2 + 2(N 2 - a)p + ftp- 1 


Biv d — model 


Phh = 0.5p _1 p|f + l-5p 3 + AHp 2 + 2{H 2 - a)p + pip- 1 



(3.53) 



5) By - models. 



By a — model 


p = pp 2 - t- 1 ^) - 7p) + t~ 2 (p - l) 2 {ap + ftp' 1 ) + Sip 


Bvb — model 


Paa = <Pp\ ~ CL^iPa - IP) + aT 2 (p - \) 2 (ap + /J^ 1 ) + Sip 


Bye — model 


Pnn = cPPn - N~Hpn - TP) + N ~ 2 (P - 1) 2 ("P + PP' 1 ) + H 


Byo — model 


Phh = <Pp 2 h - H^ipH - IP) + H~ 2 (p - l) 2 (ap + Pp- 1 ) + Sip 



(3.54) 



Here <p = Q.bp' 1 + (p- l)" 1 , ip=p(p+l)(p- l)" 1 . 
6) By i - models. 



By i a — model 


p = ^{t)p 2 mp p{ U% t] [a + $ + cm 


By ib — model 




Byic — model 


pnn = v(n)p 2 n w)pn - 'fcfff^fc + ^+ cm 


By id — model 


phh = v(h)p 2 h aH)p H + ^ + am 



(3.55) 



Here <p{z) = 0.5b- 1 + (p - I)" 1 + (p - a)" 1 ], = I*" 1 + (* - l)" 1 + (p - a)" 1 ], C(*) 

[ 7 (z-l)(p-l)- 2 + ,5z(z-l)(p-z)- 2 ]. 

3.2 Nonintegrable 0.5 - models 

Now we give some examples nonintegrable 0.5 - models. 

3.2.1 p - models 



1) A cosmology. 

i) Example 1: 

ii) Example 2: 

iii) Example 3: 

iv) Example 4: 



p = 0.5A- 1.5p 2 . 
Pad = 0.5A- 1.5/9 2 . 
Phh = 0.5A- 1.5p|. 



(3.56) 
(3.57) 
(3.58) 

N . (3.59) 

2) Pinney cosmology. It induced by the Pinncy equation. Let us present 4 examples of such 
models, 
i) Example 1: 



Pnn = 0.5A - 1.5p 2 



p = £(f)p+i- 



(3.60) 



ii) Example 2: 



iii) Example 3: 



Paa = £(a)p + — . 



phh = a n )p- 



(3.61) 



(3.62) 



iv) Example 4: 



pNN = aN)p+^- 



3) Schrodinger cosmology. For this model we give 4 submodels, 
i) Example 1: 

p = u{t)p + kp, 



(3.63) 



(3.64) 



ii) Example 2: 

paa = u(a)p + kp. (3.65) 

iii) Example 3: 

Phh = u{H)p + kp. (3.66) 
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iv) Example 4: 

Pnn = u(N)p + kp. (3.67) 

4) Hypergeometric cosmology. Let us we present 4 examples. 

i) Example 1: 

p = t-\l - + P + !)* - l]p + <*M- (3-68) 

ii) Example 2: 

Paa = a- x (l - a)- 1 ^ + p + l) a - 7 ]p a + a/3p}. (3.69) 

iii) Example 3: 

Phh =H- 1 (l-H)- 1 {[(a + l3 + l)H- 1 ]p H + al3p}. (3.70) 

iv) Example 4: 

Pnn = iV _1 (l - A0 _1 {[(a + £ + 1)W - 7W + (3.71) 

3.2.2 p - models 

1) A cosmology. 

1) Example 1: 

p = 0.5A- 1.5p 2 . (3.72) 

ii) Example 2: 

p aa = 0.5A-1.5p2. (3.73) 

iii) Example 3: 

Phh = 0.5A - l.5p 2 H . (3.74) 

iv) Example 4: 

PjV jv = 0.5A-1.5^. (3.75) 

2) Pinney cosmology. It induced by the Pinncy equation. Let us present 4 examples of such 
models. 

i) Example l: 



where £ = £(t), k — const. 
ii) Example 2: 



iii) Example 3: 



P = »+|p (3-76) 



Paa=C(a)P+^- (3-77) 



VHH=i{H)p+^. (3.78) 



iv) Example 4: 



Paw =£(A0p+^. (3.79) 

3) Schrodinger cosmology. For this model also write 5 submodels. 

i) Example l: 

p = u(t)p + kp, (3.80) 

where u = u{t), k = const. 

ii) Example 2: 

Paa = u(a)p + kp. (3.81) 
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iii) Example 3: 

Phh =u(H)p + kp. (3.82) 

iv) Example 4: 

PNN = u(N)p + k P . (3.83) 

4) Hypergeometric cosmology. Let us we present 4 examples, 
i) Example 1: 

p = t-\l - + 13 + l)t - j}p + a/3p}. (3.84) 



ii) Example 2: 

iii) Example 3: 

iv) Example 4: 



Paa = a-^l - a)- 1 {[{a + /3 + l)a- 7 ]p a + a/3p}. (3.85) 
p HH = H-\l - H)- 1 ^ + (3 + l)H- j}p H + (3-86) 
Pnn = N~\l - N)- 1 ^ + + l)N- j]p N + af3p}. (3.87) 



4 Differential equations of state: 1.0 - models 

Here we present some examples integrable and nonintegrable 1.0 - models that means models with 
the full EoS (full =1.0). 

4.1 Integrable 1.0 - models 

Here the list integrable 1.0 - models. 



Ki — model 


p pp = 6p 2 + p ' 


Ku — model 


p pp = 2p 3 + pp + a i 


Km — model 


p PP = y P - i p(pp-ap 2 -i3)+iP 3 +^ i 


Kiv — model 


Ppp = ^p 2 p + l-5p 3 + ^pp 2 + 2(p 2 -a)p+l ; | 


Ky — model 


Ppp = + ih)p 2 P - 1 (pp- ip) + p- 2 (p i) 2 («p + fo- 1 ) + ^r 1 1 


Kyi — model 


Ppp = <P(P)P 2 P - i{p)p P + P~ 2 {P - ±)~ 2 P(P - 1)(P - p)[a + Ppp' 2 + C(p)} | 



(4.1) 



Here tp(z) = 0.5b" 1 + (p - l)" 1 + (p - z)" 1 ], £(z) = [z' 1 + (z - l)" 1 + (p - z)" 1 ], C(*) - 
[7(2 — l)(p — 1) 2 + (5z(z — l)(p — z) 2 ]. All above presented Kj - models admit exact solutions. 
Now let's present some of these solutions. 

4.1.1 K/ - model 

Ppp = 6p 2 + p. (4.2) 
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4.1.2 K// - model 

Ppp = 2p 3 + pp + a. (4.3) 
This models has the following particular solutions 

V = p(p;1.5)=Tp-(2ip 2 + p)- 1 , (4.4) 

p = p(p;l) = -~ (4.5) 
P 



P = P(P;2) = ---|^7, (4.6) 
p p 3 + 4 



p = ^ ;3) = 7TI' ^ + 2o P 3- 8 o ' ( 4J ) 

1 6p 2 (p 3 + 10) 9p 5 (p 3 +40) 

P ~ PW J p p 6 + 20p 3 -80 p 9 +60p 6 + 11200' 1 j 

p = p(p; 0.5e) = -e^> (4.9) 

and so on. Here 

1> = (In 4>) p , 4>{p) = C 1 A J (-2- 1 / 3 p) + C 2 Bi(-2-^ 3 p), (4.10) 
where Cj = consts and Bi(x) are Airy functions. 



Ppp = -Pp - ^(Pp - «P 2 - /3) + 1P Z + -• (4.H) 



4.1.3 K in - model 

-Pp - ~(Pp - ap 2 -P) + 1P Z + -• 
P p P P 

This equation admits the infinite number exact solutions. For example, it has the following par- 
ticular solutions (see e.g. [I]) 

p = p(p;p, -pre 2 , A, -Are 4 ) = re, (4.12) 
P = p(p;0,-p,0,pk) = «p, (4.13) 

p = p(p; 2re + 3, -2re+ 1,1,-1)= ^ + * , (4.14) 

P + re + 1 

p = p(p;p, 0,0, -pre 3 ) = ntfp, (4.15) 

p ee p(p;0, -2re,0,4rep- A 2 ) = p[re(lnp) 2 + A In p + p], (4.16) 

p = p(p;-^ 2 A,0^ 2 (A 2 -4rep),0)= ^— , (4.17) 

rep^ + Ap y + p 

p = p(p;0.5e) = -e 1 (ln ( /») p (4.18) 

and so on. Here 

<t>{p)=p v [C l J v {Q + C 2 Y v {C,)] ) (4.19) 
where Cj = consts, £ = ^/e^p, f = 0.5aei and ^(C)>^(C) are Bessel functions. 

4.1.4 K/y - model 

Ppp = :r p p + L5p3 + 4 ^ 2 + 2 (-° 2 - a ^P + ( 4 ' 2 °) 

zp p 

This equation has the following particular solutions (see e.g. [4]): 

p = p(p;±2,-2)=±p- 1 , (4.21) 

p ee p(p;0,-2) = -2p, (4.22) 

p ee p(p; 0, -3~ 2 2) = -3 _1 2p, (4.23) 

p ee p(p;-m,-2(m-l) 2 ) = -[lnff m _i(p)] p , (4.24) 

p ee p(p) = 2iTT- a - 5 e p2 [iC + erfc{ip)]-\ (4.25) 

p ee p(p)=2n-°- 5 e-P 2 [C-erfc(p)}-\ (4.26) 

p ee p(p; 0.5e) = - £l (In <(>)p (4.27) 
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and so on. Here 

<j>{p) = [C\U(C,2°- 5 p) + C 2 l/(C,2°- 5 p)]e a5£p2 , (4.28) 

where Cj = consts, C, = a + 0.5e, ^ = 0.5aei and U, V) are parabolic cylinder functions, H m 
are Hermite polynomials, erfc is the complementary error function (for detail see e.g. [1]). 

4.1.5 Ky - model 

Ppp = + -^-r)^ - -(P, - TP) + P~ 2 (p - 1) 2 («P + /fr -1 ) + 6p(p+1) . (4.29) 
2p p — 1 ^ p p — 1 

This equation admits the following particular solutions (see e.g. [4]): 

(4.30) 
(4.31) 
(4.32) 
(4.33) 
(4.34) 
(4.35) 



p 


= p{p; 0.5, -0.5p 2 , k(2 - p), -0.5k 2 ) = up - 


\- p, 


p 


= p(p; 0.5, K 2 p, 2np, p) = n(p + k) _1 , 




p 


= p(p; 0.125, -0.125, -np, p) = (/c + p)(k - 


pT 


p 


= p(p;p, -0.125, -pK 2 ,0) = 1 + np - 5 , 




p 


= p{p; 0, 0, p, -0.5p 2 ) = Ke w , 




p 


= P(fi) = -eipQn<l>)p 





and so on. Here 



<t>{p) = {e 2 p)- v [C x M R ^{e 2 p) + C 2 W K Je 2 p)}e - 5t2P , (4.36) 
where Cj = consts and M, W are Whittaker functions (for more exact details see e.g. [3]). 



4.1.6 Kur - model 



1 



p = 0.5 I - + — - + — — \p 2 -[- + — — + — J Pp 

\p p-1 p-p) p \p p-1 p-pj 

+ p- 2 (p - l)- 2 p( P - l)(p - p) [a + /?PP~ 2 + 7(P - 1)(P - I)" 2 + MP - 1)(P - P)~ 2 ] • (4-37) 
This equation is integrable and has the following particular solutions (see e.g. [3]): 



p 


= p(p; p, — pK 2 , 0.5, 0.5 — p(« — l) 2 ) = rep, 


(4.38) 


p 


= p(p;0,0,2,0) = K p 2 , 


(4.39) 


p 


= p(p; 0,0, 0.5, -1.5) = up- 1 , 


(4.40) 


I' 


= p(p;0,0,2,-4) = «p- 2 , 


(4.41) 


p 


= p(p; 0.5(k + p) 2 , -0.5, 0.5(p - l) 2 , 0.5k(2 - k)) = p(re + pp) -1 , 


(4.42) 


p 


= p(p; 0.5k 2 , -0.5k 2 , 0.5p 2 , 0.5(1 - p 2 )) = p 5 


(4.43) 



and so on (for more exact details see e.g. [4]). 

4.2 Nonintegrable 1.0 - models 

1) Example 1: 



p w = 0.5A-1.5p 2 . (4.44) 

Ppp = £(p)p+4' (4-45) 

P PP = u (p)p + k P, (4-46) 

Ppp = p-\l - Pallia + /3 + l)p- 7 ]p + (4.47) 



2) Example 2: 

3) Example 3: 

4) Example 4: 

5) Example 5: 

Ppp = n{n - 1) V^V^ 2 - (4.48) 

This model has following solution 

p = -Ap n . (4.49) 
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5 D - models 



In this section we would like to consider the so-called D - models, where D = D(p, p) is an arbitrary 
function of p and p. Here some examples of such models. 

5.1 Integrable D - models 
1) Di - models. 



Dia — model 


D = 6D 2 + t | 


Djb — model 


D aa = 6D 2 + a | 


Die — model 


D NN = 6D 2 + N | 


Did — model 


D HH = 6D 2 + H | 


Die — model 


D pp = 6D 2 + p | 


Dip — model 


D pp = 6D 2 +p j 



(5.1) 



2) Dn - models. 



Diia — model 


D = 2D 3 + W + a j 


Dub — model 


D aa = 2D 3 + aD + a | 


Due ~ model 


D NN = 2D 3 + ND + a | 


Dn d — model 


Dun = 2D 3 + HD + a | 


Due — model 


D pp = 2D 3 + pD + a | 


Due — model 


D pp = 2D 3 +pD + a | 



(5.2) 



3) Dm - models. 



Dm a — model 


D = D- 1 ^ 2 -t-\D - aD 2 - /3) +jD 3 + SD- 1 \ 


Dihb — model 


D aa = D- X D\ - a-^Da - aD 2 - 0) + 1 D 3 + SD' 1 \ 


Dine — model 


D NN = D^Djf - N-^Dn - aD 2 - 0) + jD 3 + SD' 1 \ 


Dihd — model 


Dhh = D~ 1 D 2 H - H-\D H - aD 2 - p) + 1 D 3 + SD' 1 \ 


Dihe — model 


D pp = D- x D 2 p - p-\D p - aD 2 - /?) + 7 D 3 + SD' 1 \ 


Dihf — model 


D pp = D- x D 2 p - p~\D p - aD 2 - (3) + 1 D 3 + SD' 1 | 



(5.3) 



4) Div - models. 
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Diva — model 


D = 0.5D- 1 !) 2 + 1.5£> 3 + AtD 2 + 2{t 2 - a)D + fiD^ 1 | 


Divb — model 


D aa = {).5D- l D 2 a + 1.5D 3 + AaD 2 + 2(a 2 - a)D + (3D- 1 \ 


Dive — model 


D NN = O.SD^Djj + 1.5D 3 + AND 2 + 2(N 2 - a)D + (3D- 1 \ 


Divd — model 


D HH = Q.hD^D 2 ] + 1.5D 3 + AHD 2 + 2{H 2 - a)D + (3D- 1 | 


Dive — model 


D pp = Q.hD^Dl + 1.5D 3 + ApD 2 + 2(p 2 - a)D + (3D' 1 \ 


Divf — model 


D pp = Q.hD^Dl + 1.5D 3 + ApD 2 + 2{p 2 - a)D + (3D' 1 | 



5) Dv - models. 



Dva — model 


D = <t)b 2 - t-\D - 7 L>) + t~ 2 (D - l) 2 (aD + (3D- 1 ) + Sip \ 


Dvb — model 


D aa = cf>D 2 a - a-^Da - jD) + a- 2 (D - l) 2 (aD + /3D- 1 ) + 5^ \ 


Dvc — model 


D NN = 4>D 2 N - N~\D N - jD) + N- 2 (D - l) 2 (aD + (3D- 1 ) + 5i/> \ 


Dvd — model 


Dhh = 4>D 2 H - H-\D H - lD) + H- 2 (D - \) 2 {aD + (3D- 1 ) + 5^ \ 


Dve — model 


D pp = <f>D 2 p - p-\D p - 7 £>) + p- 2 (D - l) 2 (aD + (3D- 1 ) + <ty \ 


Dv f — model 


D pp = <^D 2 V - p~\D p - 1 D)+p- 2 {D - l) 2 (aD + (3D- 1 ) + 5^ \ 



Here = 0.5D- 1 + (D - l)" 1 , ip = D(D + 1)(D - l)" 1 . 



6) Dvi - models. 



Dvia — model 


b = V (t)b 2 mb yy'ia + n + at)} | 


Dvib — model 


D aa = v(a)D 2 a - i(a)D a ^gg^ [a + % + ((a)} | 


Dvic — model 


D NN = cp(N)D 2 N - £(N)D N - D %%% N) [* + §f + ((N)} | 


Dvib — model 


Dun = <P{H)D% S{H)D H D %^% H) [a + §f + ((H)] \ 


Dvie — model 


D pp = ^p)D 2 p i{p)D p D{ X~ii P) [ Q + W + C(P)] | 


Dvif — model 


D pp = ^(p)D 2 p - ttp)D p D( X-^ P) i a +W + C(P)] | 



Here 

<p{z) - Q.51D- 1 + (D - l)- 1 + (D - z)- 1 }, 

£{z) = z-' + iz-iy' + iD-z)- 1 , 

C(z) - 7 (z-l)(£>-l)- 2 + 5z(z-l)(D-z)- 2 . 
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5.2 Nonintegrable D - models 

Nonintegrablc D - models can be induced by some ODE's which are nonintegrable e.g. as in our 
previous papers [I]-j2]- 

5.3 Solutions 

It is important that all above presented D - models admit exact solutions. Let's consider here one 
example. Let D(p,p) has the form 

D = 0.5(p + p). (5.10) 
Consider the Dvid - model [see the system (5.6)]. This model in fact has the form 

p = -2H-3H 2 , (5.11) 
p = 3H 2 , (5.12) 

D HH = V {H)D% - C(H)D H - D{D H ^H-lf H) t Q + if + C(#)] (5-13) 
p = -3H(p + p), (5.14) 

Our aim is solve this system. For the particular case (5.10), the equation (5.13) has the following 
particular solution [4] 

D = kH 2 , (k = const). (5.15) 



From (5.11)-(5.12) we get 
which has the solution 



H = -kH 2 . (5.16) 



H = -. (5.17) 



6 Artificial two-dimensional models 

FRW cosmological models that we considered above are one-dimensional. But for some reason we 
would like to have two-dimensional cosmological models in FRW spacetime. And there are (may 
be) no legal ways to construct two-dimensional models starting from one-dimensional models. If 
so let's try to use the "illegal" ways e.g. introducing the artificial "coordinate". As such artificial 
coordinate we can use e.g. one of physical parameters of the original model e.g. the cosmological 
constant (A). So we have may be two coordinates: one legal coordinate - t (time) and one "illegal" 
coordinate - A. Now we are ready to write our artificial two-dimensional cosmological models in 
FRW spacetime. 

1) As an example, consider the following model e.g. for the scale factor [a a — da/dA, a t — 
da/dt etc]: 

at = 0.75(a 2 ) A (6.1) 

or its twin 

a A = 0.75(a 2 ) t (6.2) 

which arc known to develop shocks [Exercise 1: What means these shocks for the dynamics of the 
universe?]. Eqs. (6.1) and (6.2) are nothing but the dispersionless Korteweg-de Vries equations 
(dKdVE) or Riemann equations. It is well-known that the equations (6.1) and (6.2) have the 
following solutions 

a(A,t) = /i(A-0.75ai) (6.3) 

and 

a(t,A) = h(t - 0.75aA), (6.4) 

respectively. Here h is an arbitrary function. Also we see e.g. from the solution (6.3) that 
the velocity of a point of the wave, with constant amplitude a, is proportional to its amplitude 
leading to the "breaking" of the wave. Also we note that the wave also develops discontinuities in 
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its evolution [Exercise 2: What means these discontinuities for the dynamics (evolution) of the 
universe?]. Let's now we give the following particular solutions of the equations (6.1) and (6.2): 

a=(/3i+/3 2 A)(/3 3 -1.5/3 2 *)- 1 (6.5) 

and 

a = (/?i+/3 2 t)(/33-l-5/3 2 A)- 1 , (6.6) 

respectively, where = consts. 

2) Our second example is the Euler-Tricomi equation 

att = ta AA . (6.7) 
It is known that this equation has the following particular solution 

a = a(3A 2 + i 3 ) + /3(A 3 + t 3 A) + 5{6tA 2 + < 4 ). (6.8) 
The twin of the equation (6.7) reads as 

a AA = Aa tt - (6.9) 

3) Let's now present the Dym equation 

a t = {a-°- 5 ) AAA (6.10) 

and its twin 

a A = (a-°- 5 ) ttt . (6.11) 

4) Our next example is given by [5] 

(lno) t i = OAA (6-12) 

and its twin 

(lna) A A = a tt . (6.13) 

5) Our last example is given by 

a t = fJ-a AA . (6-14) 

It is the heat equation and it is well-known that it has the following fundamental solution [as 
a(A,t = 0) = 5(A)] 

a-(47r^)-°- 5 e -°- 25 ^ lt " A2 . (6.15) 



Note that the twin of the equation (6.14) looks like 

a A = fxa tt - (6.16) 

7 Conclusion 

It is important to study the physical and mathematical nature of FRW cosmological models as they 
play a crucial role in modern cosmology. In this work, a new classes integrable and nonintegrable 
FRW cosmological models were proposed. To construct integrable models we implant integrable 
equations, in our Painleve equations, into the body of the original gravitational equations. It is 
remarkable that all proposed integrable models admit exact solutions. For some of them, exact 
solutions are presented. Finally in Appendix, we extend the obtained results for g-essence models 
and its two reductions: k-essence and f-essence. 
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8 Appendix. Integrable and nonintegrable g-essence mod- 
els and their k-essence and f-essence reductions 

For g-essence we have 

p = K, p = 2XK x +YK Y -K, (8.1) 

where K is the Lagrangian for g-essence, X and Y are the kinetic terms for the scalar and spinor 
fields, respectively. Below we give some examples of integrable models for g-essence. 



2) Ej - models. 



Ei a — model 


K = 6K 2 + t 


Ejb — model 


K aa = 6K 2 + a | 


Ejc — model 


K NN = 6K 2 + N | 


Ej£> — model 


K HH = 6K 2 + H\ 


Eje — model 


K xx = GK 2 + X | 


Eif — model 


Kyy = 6K 2 + Y | 



(8.2) 



2) Ejj - models. 



En a — model 


K = 2K 3 + tK + a | 


Ei ib — model 


K aa = 2K 3 + aK + a \ 


Enc — model 


K NN — 2K 3 + NK + a | 


En D — model 


K HH = 2K 3 + HK + a \ 


En E — model 


K XX = 2K 3 + XK + a\ 


Ei if — model 


K YY = 2K 3 + YK + a | 



(8.3) 



3) Em - models. 



EniA — model 


K = K- 1 k 2 -t- 1 (k-aK 2 -f3)+- f K 3 + 5K- 1 | 


Ei 1 1 b — model 


K aa = K- x K 2 a - a-\K a - aK 2 - /3) + jK 3 + SK' 1 \ 


Einc — model 


K NN = K^K 2 N - N- X {K N - aK 2 - f3) + jK 3 + SK' 1 \ 


En id ~ model 


K HH = K- X K 2 H - H-HKh - aK 2 - f3) + jK 3 + SK' 1 \ 


Eiiie — model 


Kxx = K- X K 2 X - X- X {K X - aK 2 - p) + ~{K 3 + 5K- 1 \ 


Eiiif — model 


K YY = K^K 2 - Y~ 1 (K Y - aK 2 - j3) + jK 3 + 5K' 1 \ 



(8.4) 
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4) Ej V - models. 



Eiva — model 


K = O.BK^K 2 + 1.5K 3 + UK 2 + 2(t 2 - a)K + /3R- 1 \ 


Eiv B — model 


K aa = 0.5K~ 1 Kl + 1.5K 3 + 4aK 2 + 2(a 2 - a)K + jiK' 1 \ 


Ejvc — model 


K NN = 0.5K~ 1 K% + 1.5K 3 + ANK 2 + 2(N 2 - a)K + pK' 1 | 


Ejv d — model 


Khh = Q.bK^Kjj + 1.5K 3 + 4HK 2 + 2{H 2 - a)K + /3R- 1 \ 


Eiv e — model 


K X x = ^.hK- x K\ + 1.5K 3 + 4XK 2 + 2(X 2 - a)K + /3R- 1 \ 


Eiv d — model 


K YY = O-ZK-iK 2 + 1.5K 3 + 4YK 2 + 2(Y 2 - a)K + PR- 1 \ 



5) Ev - models. 



Eva — model 


K = 4>K 2 - t-^K - jK) + t- 2 (K - l) 2 (aK + PR- 1 ) + 5%j} \ 


Ev b — model 


K aa = 4>K 2 a - a-\K a - jK) + a- 2 (K - l) 2 (aK + PK' 1 ) + 8^ \ 


Eve — model 


K NN = cj)K 2 N - N-^Kn - <yK) + N- 2 (K - l) 2 (aK + /3R- 1 ) + 5^ \ 


Ev d — model 


K HH = 4>K 2 H - H~\K H - jK) + R- 2 {K - \) 2 {aK + PR- 1 ) + 5^ \ 


Ev e — model 


K X x = 4>K 2 X - X- X {K X - jK) + X- 2 (K - l) 2 (aK + PR- 1 ) + Sif) \ 


Ev d — model 


K YY = 4>K 2 - Y- X (K Y - jK) + Y- 2 (K - \) 2 {aK + PR- 1 ) + Sip | 



Here = 0.5D- 1 + (R - l)" 1 , V = K(K + 1)(K - l)" 1 . 



6) Evi - models. 



Evia — model 


K = V (t)K 2 t(t)K K{ ^l[«r t] [a + || + C(t)] | 


Evi b — model 


K aa = v(a)K 2 t(a)K a - K %^ a) [^ + ^ + C(«)] | 


Evic — model 


K NN = V {N)K 2 N d{N)K N - K{ ^% N) {a + ^ + Q{N)] \ 


Evi d — model 


K HH = f(H)K 2 H i{H)K H - K{ ^% H) [a + ^ + ((H)] \ 


Evie — model 


K XX = <p(X)K% i(X)K x ~ + f| + CPO] | 


Evif — model 


K YY = V(Y)K 2 i(Y)K Y K( yhy% Y) [" + ^ + COO] | 
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